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Abstract 

Effects of interlayer coupling and spatial anisotropy on spin-wave excitation spectra of a three- 
dimensional spatially anisotropic, frustrated spin-^ Heisenberg antiferromagnet (HAFM) is inves- 
tigated for the two ordered phases using second-order spin-wave expansion. We show that the 
second-order corrections to the spin-wave energies are significant and find that the energy spectra 
of the three-dimensional HAFM shares similar qualitative features with the energy spectra of the 
two-dimensional HAFM on a square lattice. We also discuss the features that can provide exper- 
imental measures for the strength of the interlayer coupling, spatial anisotropy parameter, and 
magnetic frustration. 
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The intriguing properties of many layered antiferromagnets and recently discovered iron 
based superconductors have created considerable interest in two and three dimensional spin-^ 
Heisenberg antiferromagnet (HAFM) with competing interactions.^ For the last two decades 
the properties of quantum spin-^ Heisenberg antiferromagnet (HAFM) with nearest neighbor 
(NN) Ji and next nearest neighbor exchange interactions (NNN) J2 on a square lattice have 
been studied extensively by various analytical and numerical techniques.-i^ Earlier studies 
on the Ji — J2 model have shown that the ground state of the HAFM is antiferromagnet- 
ically (AF) ordered at low temperatures. Addition of J2 interactions destabilize the AF 
order and for a critical value of the frustration parameter rj = J2/ Ji a. quantum disordered 
phase emerges. With increasing values of J2 there is a second quantum phase transition 
from the disordered phase to a columnar antiferromagnetic (CAF) stripe phase. Using 
experimental techniques such as nuclear magnetic resonance, magnetization, specific heat, 
and muon spin rotation measurements, properties of layered magnetic materials Li2VOSi04, 
Li2VOGe04, VOM0O4, and BaCdVO(P04)2 have been studied." These studies have shown 
that these compounds have significant couplings between NN and NNN neighbors. More- 
over, for Li2VOSi04, a layered material that can be described by a square lattice Ji — J2 
model with large J2 the interlayer coupling J±/Ji ~ 0.07 is not negligible.^ Due to a finite 
interlayer magnetic coupling J±, these experimental systems are quasi 2D. 

Experimentally one of the most direct ways to probe the magnetic excitation spectra is 
with inelastic neutron scattering (INS). For example, INS experiments on S* = 1/2 2D anti- 
ferromagnets Sr2Cu304Cl2, Cu(DCOO)2-4D20 (CFTD), and La2Cu04 (LCO) have revealed 
that the magnon energy at the wave- vector (tt, 0) is depressed by 7% relative to the energy 
at the wave- vector (7r/2, 7r/2).^°"-^ Another example is the recently discovered iron pnictide 
superconductors.^'^ In the parent phases of these materials magnetic excitations have shown 
to play an important role in the superconducting state.—"— INS measurements on CaFe2As2 
and SrFe2As2 have revealed a very large difference between Ji and J{.— ~— 

Theoretically the spin-wave spectra and the low-temperature magnetic phase diagram 
have been obtained for the spatially anisotropic Ji — J[ — J2 HAFM model on a square 
lattice with NN exchanges Ji along the x axis, J[ along the y axis, and NNN interactions 
J2 along the diagonals in the xy plane.^ Recent experiments on iron-based superconduc- 
tors such as undoped iron oxypnictides reveal that the electronic couplings are more three 
dimensional than in the cuprate superconductors.—"— With decrease in temperature most 
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undoped iron-pnictide superconductors show a structural transition from a tetragonal para- 
magnetic phase to a orthorhombic phase. In the 122 materials a three dimensional (3D) 
long-range antiferromagnetic order develops simultaneously. The magnetic phase diagram 
of the 3D Ji — J[ — J2 — J± HAFM model with interlayer coupling has also been recently 
investigated.-i2^'2^. The present work is a sequel of Ref. l3| by the author. Here we study 
the spin-wave energy dispersion of the two ordered phases of the model taking into account 
the first (1/5*) and second-order (1/5*^) corrections to the Hamiltonian. The details of the 
model and the derivations along with many references can be found in Ref. [sj. 

The Heisenberg Hamiltonian for a spatially anisotropic, frustrated spin-| HAFM on a 
cubic lattice with four exchange AF interactions between spins: Ji along the x (row) direc- 
tion, J[ along the y (column) direction, J2 along the diagonals in the xy plane, and with 
interlayer coupling J± is described by 

£ labels the layers, i runs over Nl lattice sites, i + 6x {Sx = ±1) and i + 6y {6y = ±1) are 
the NN interactions to the i-th site along the row and the column direction. The third 
term represents the NNN interaction, which are along the diagonals in the xy plane and the 
last term is for the NN coupling between the layers. ( = J[/Ji measures the directional 
anisotropy, t] = J2/J1 is the magnetic frustration between the NN (row direction) and 
NNN spins, and 6 = J±/Ji is the interlayer coupling parameter. At zero temperature 
the classical ground states for the anisotropic model are the Neel or the AF (tt, tt, tt) state 
and the columnar antiferromagnetic (CAF) (vr, 0, vr) state (the classical transition occurs at 
Vc = C/2). 

For the AF ordered phase, NN interactions are between A and B sublattices and NNN 
interactions are between A-A and B-B sublattices. The Hamiltonian in Eq. [T]takes the form: 



This Hamiltonian is mapped to bosonic creation and annihilation operators a\a and h\h 
using Holstein-Primakoff transformations keeping only terms up to the order of l/S"^. Fur- 
thermore, the Hamiltonian is diagonalized in the momentum space using the Bogoliubov 



transformations. In powers of 1/S the Hamiltonian can be written as^ 



H = --NJ,S'z{l + 
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+ Ho + H1 + H2 + 
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^ Cik (ala^ + f3lM + Csk UIpU + /3-k«k) + 



In the above equation, 



^k = (l-7^)^/^ 

7ik = [cos(fca;) + Ccos(fcj^) + 6 cos(fc^)]/(l + C), 72k = cos(fca;) cos{ky), (7) 
7k = 7ik/ /tk, /tk = 1 



1 + '^-^^ ■ 1 + C 

The first term in Hq is the zero-point energy and the second term represents the excita- 
tion energy of the magnons within hnear spin- wave theory (LSWT). The part Hi and if 2 
correspond to 1/S and l/S"^ corrections to the Hamihonian. 

The spin-wave energy E^^ for magnon excitations, measured in units of 715*^(1 + () up 
to second order in l/S" for the AF-phase is given as^ 



Ed"" = E^ + — -Ak + 



{2S) " (25)2 
where the second-order self-energy is given by 

2 \2 
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k,p,q,[k+p-q] 



IK 
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k,p,q,[k+p-q] 



.u-Ep-E^- f;k+p-q + iO+ u + Ep + E^ + f;k+p-q + ^0- 
The coefficients ^k, -Bk, Cik are given in Ref. | 



(9) 



4 



On the other hand, the Hamiltonian describing the CAF phase is 
H = Ji ^ ■ Sf^^^ £ + -J[ ^ I^S^^ ■ ^fj^Sy/ + ■ ^%5y,e. + J2 ^ S^^ ■ S 



B 



+ J^J2^t,i-Sf,e+i- (10) 

i,e 

The quasiparticle energy E^^^, measured in units of JiSz{l + 2r]) up to second order in 1/S 
for the CAF phase is 



/2 



s2Hk,K)-^J- (11) 



k ^ 25 (25)2 
The structure factors 7^,^, in this case are 

7ik = [cos{k^){l + 2r]cos{ky)) + Scos{h)]/{l + 2r]), ^2^ = cos{ky), 
7k = 7;k/4, < = l- Y^(l - 7k) + (12) 

4 = [i-7LT/^- 

The different coefficients can be found in Ref. [sj. 

Spin-wave energy dispersions for the two ordered phases with different values of (, t], 
and 6 are obtained numerically from Eqs. |8] and [TTJ Especially to obtain the second order 
correction terms we sum up contributions from A^^^ points of p and points of q in the 
first BZ. For our results we use Nl = 18 lattice sites. 

Figured] shows the spin-wave energies E^^ / Ji for the AF phase (panels a and b) with ^ = 
1, 5 = 0.1, 7] = and for the CAF phase (panels c and d) with ^ = 0.9, 5 = 0.1, = 1 obtained 
from LSWT (solid lines), with first-order corrections (dotted lines), and with second-order 
corrections (long dashed lines). Momenta {kx,ky) are varied from (0,0) to (vr, tt) in the first 
BZ with kz fixed at vr or as (0,0, vr) — )■ (7r,0,7r) — )• (7r/2, 7r/2, tt) — )■ (7r,7r,7r) — )• (7r,0,7r) 
(panels a and c) and (0, 0, 0) -> (vr, 0, 0) -> {n/2, it/2, 0) -> (vr, vr, 0) -> (tt, 0, 0) (panels b and 
d). Panels (a) and (b) have Goldstone modes at the wave- vectors k = and k = (vr, vr, vr) 
whereas for panels (c) and (d) the Goldstone modes are at k = and k = (7r,0,7r). The 
local energy minima are at (vr, vr, vr) for the AF phase and at (vr, 0, vr) for the CAF phase. For 
both the phases the first and second order corrections significantly increase the spin-wave 
energy from that obtained with LSWT. However, for the CAF phase l/S"^ corrections lower 
the spin-wave energy from the first-order corrections. 
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FIG. 1: (Color online) Spin-wave energies / Ji for the AF phase (panels a and b) with = 
1,5 = 0.1, and r] = and for the CAF phase (panels c and d) with C = 0.9,6 = 0.1, and rj = 1 
obtained from LSWT (solid lines), with first-order corrections (dotted lines), and with second- 
order corrections (long dashed lines). Momenta {kx,ky) are varied from (0,0) to (7r,7r) in the first 
BZ keeping = tt (panels a and c) and k^ = (panels b and d). The local energy minima are 
at (vr, vr, tt) for the AF phase and at (tt, 0, tt) for the CAF phase. For both the phases the first 
and second order corrections significantly increase the spin-wave energy from the results obtained 
using LSWT. For the CAF phase l/S"^ corrections lower the spin- wave energy from the first-order 
corrections. 

The spin-wave energies for AF phase E^^ / Ji is plotted in Fig. [2] for different values 
of C,1],6 in the BZ. The dispersion along the (7r/2, 7r/2, tt) — (7r,0,7r) and (7r/2, 7r/2, 0) — 
(71,0,0) are flat within LSWT and 1/S correction (see Fig. [1] for example). The second 
order corrections make the magnon energies at (tt, 0, tt) or (tt, 0, 0) smaller than the energies 
at (vr/2, 7r/2, vr) or (vr/2, 71/2, 0). Also we find that with increase in NNN frustration the dip 
in the spin-wave energies at (vr, 0, vr) and (vr, 0, 0) increases. This can provide an experimental 
measure of the strength of the NNN frustration. These features are qualitatively similar to 
the recently obtained spin-wave energy dispersion for the frustrated HAFM on the spatially 
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FIG. 2: (Color online) Plot of spin-wave energy dispersion (with second-order corrections) / Ji 
for the AF phase with = 1, 6 = 0.1 (panels a and b), 0.3 (panels c and d), rj = 0, and 0.2. 
kz = IT for panels a and c and kz = for panels b and d. Minimum of energy occurs at (vr, vr, vr) for 
the AF phase (the narrow E^^ = window near (vr, 7r,-7r) is due to numerical inaccuracy). With 
increase in NNN frustration the dip in the spin-wave energies at (vr, 0, vr) and (vr, 0, 0) increases. 
This can provide an experimental measure of the strength of NNN frustration. The strength of the 
interlayer coupling 6 an be measured from the spin- wave energy at (0,0, vr) (see text). 

anisotropic square lattice .- 

Figure E] displays the energy dispersion for the CAF phase with ( = 0.9, S = 0.1 (panels 
a and b), 0.3 (panels c and d), rj = 1 and 0.6. We find three energy peaks, the maxima 
being at (vr/2,0,vr) and (vr/2,0,0). With decrease in frustration the system approaches the 
phase transition point, so the spin- wave energy diminishes. 

For 5 = 0, we have a Goldstone mode at the wave-vector k = (0, 0, vr) for both the 
AF and CAF phases. But presence of a finite interlayer coupling 6 opens up a gap at 
these wave- vectors (see Figs. [2] and [2] (panels a and c)). This gap increases with increase 
in 6. We find that the spin-wave energy E^^ / Ji for the AF-phase obtained from LSWT 
(using Eq. [8]) at the wave-vector (0, 0, vr) in the BZ is 2^J5{1 + C)/[l + 5/(1 + C)], which is 




FIG. 3: (Color online) Spin-wave energy dispersion (with second-order corrections) / Ji for 

the CAF phase with ( = 0.9, r] = 1,0.6, and 5 = 0.1 (panels a and b), 0.3 (panels c and d). 
Minimum of energy occurs at (vr, 0, vr) for the CAF phase. The energy gap at the wave-vector 
(0, 0, it) increases with increase in 6 (see panels a and c) 

independent of frustration t] (see Fig. [2|i). As an example, within LSWT ioi ( = 1 we find 
E^^(0,0,7r)/Ji ^ 0.85 with 6 = 0.1 (Fig. On the other hand, for the CAF phase we 
find the energy ^^^^(0, 0, 7r)/Ji (from LSWT) to be 2^5(1 + 2r])/[l + 5/(1 + 2r])], which 
is independent of the anisotropy parameter ( (see Fig. Ht). In this case we find for rj = 1, 
£]CAF^Q^ 0, n)/ Ji ~ 1.1 with 6 = 0.1 (see Fig.[Th). However, first and second order corrections 
slightly increase these LSWT predictions. The strength of the interlayer coupling 6 can be 
measured experimentally from the energy at this k = (0, 0, vr) value. 

In Fig. Hlwe show the effects of spatial anisotropy on the two ordered phases of the model. 
For the AF phase with increase in directional anisotropy (i.e., with lowering values of () the 
spin-wave energy decreases. The dip in the energy at (vr, 0, 0) or (tt, 0, vr) can provide an 
experimental measure of For the CAF phase the trend is opposite, i.e. E^^^ / Ji increases 
as C is lowered from 1 to 0.6. From panel c we find that at the wave-vector (0,0, vr), 
£''^^^(0, 0, 7r)/Ji is independent of C- This is similar to i?'^^(0, 0, 7r)/Ji being independent 
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FIG. 4: (Color online) Effects of spatial anisotropy C on the energy dispersion for the AF and 
the CAF phases are shown. For the AF phase with increase in directional anisotropy (i.e. with 
lowering values of () the spin-wave energy decreases. The dip in the energy at (vr, 0, 0) or (vr, 0, vr) 
can provide an experimental measure of (. For the CAF phase the trend is opposite, i.e. E^^^ / J\ 
increases as Q is lowered from 1 to 0.6. £^^^^(0, 0, 7r)/Ji is independent of C at the wave- vector 
(0, 0, vr) (panel c). 

of 7] (see panels a and c of Figl2]) 

The present work is a sequel of an earlier theoretical paper where the author studied the 
magnetic phase diagram of a spatially anisotropic, frustrated spin-i HAFM on a 3D lattice 
with interlayer coupling.^ In the present work, we have studied the effects of interlayer 
coupling and spatial anisotropy on the spin-wave energies of the two long-range ordered 
phases (antiferromagnetic Neel and antiferromagnetic columnar stripe). We have shown 
that the second-order corrections to the energy spectra are significant. Also we have found 
that our obtained energy spectra for the three-dimensional HAFM model (with interlayer 
coupling) shares similar qualitative features with the energy spectra for the two-dimensional 
square lattice. Finally, we have provided a few key features in the energy spectra that can be 
measured experimentally (e.g., with neutron scattering experiments). These measurements 
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can provide us information on the strength of the interlayer couphng, spatial anisotropy 
parameter, and magnetic frustration. 
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